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ION AND ELECTRON TRANSPORT IN AN 
NSTAR-DERIVATIVE ION THRUSTER 


John E. Foster 

National Aeronautics and Space Administration 
Glenn Research Center 
Cleveland, Ohio 44135 


SUMMARY 

Diffusion of electrons and ions to anode surfaces between the magnetic cusps of an NASA Solar 
Electric Propulsion Technology Application Readiness ion thruster has been characterized. Ion flux 
measurements were made at the anode and at the screen grid electrode. The measurements indicated that 
the average ion current density at the anode and at the screen grid were approximately equal. Additionally, 
it was found that the electron flux to the anode between cusps is best described by the classical cross-field 
diffusion coefficient. 


NOMENCLATURE 


B transverse magnetic field 

d eff constant of proportionality 

D e electron diffusion coefficient, B=0 

D ± electron cross-field diffusion coefficient 

e elementary charge of an electron 

E s sheath electric field 

J e electron cross-field current density 

J scm grd ion current density at the screen grid 

J , o ion current density at wall probe B 

k Boltzmann's constant 

m mass of an electron 

ll e electron number density 

r radial displacement 

T A electron temperature at probe A 

T B electron temperature at probe B 

T e electron temperature 

A V voltage drop across plasma sheath 

fi electron mobility 

V e total electron collision frequency 

W electron cyclotron frequency 
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INTRODUCTION 


Controlling ion and electron transport to discharge chamber surfaces is a key aspect of ion thruster 
design. 1,2 Indeed, at a given thruster power, performance can be optimized by maximizing the fraction of 
discharge ions that reach the ion extraction optics. Maximizing the extracted ion fraction can be achieved 
by minimizing ion losses to the anode. Reducing the flow of ions to surfaces other than the extraction 
optics using a magnetic field is usually not practical due to the very strong field strengths required to 
magnetize ion motion. Instead, thruster magnetic circuits, such as the ring cusp, utilize magnetic field 
geometries that contain electrons in a manner such that the discharge produced is located just upstream of 
the ion optics. 3 The close proximity of discharge to the collection optics increases the likelihood that an ion 
generated in this volume will be extracted. 

It has been suggested that many ion thruster discharge chamber configurations operate with 
preferential ion flow to the ion optics. 1 Whether or not the ring cusp magnetic circuit of the NASA Solar 
Electric Propulsion Technology Application Readiness (NSTAR) 4 ion thruster investigated in this work can 
effect ion motion to the anode and thereby improve ion flow to the optics is a question that is addressed in 
this work. An understanding of the nature of ion flow from the discharge to the ion optics can be used as a 
guide in both optimizing the performance of this thruster and also increasing the thruster's power range. 

Unlike ion motion, electron motion is severely affected by the presence of the inter-cusp magnetic 
field. In this work, electron diffusion to the anode surfaces between the magnetic cusps was also studied. 
The magnetic field between cusps provides a fairly high impedance path for electron flow to the anode, 
thereby reducing electron losses to the anode. The character of the cross-field electron diffusion coefficient 
is a measure of this impedance. A large cross-field diffusion coefficient would lead to a lower discharge 
voltage and consequently, a reduction in the bulk plasma ionization rate. Previous investigations of electron 
diffusion between cusps have assumed anomalous Bohm diffusion/ In this work, the nature of the cross- 
field electron diffusion coefficient in a NSTAR-derivative ion thruster is determined. 


EXPERIMENTAL APPARATUS 

The ion thruster used in this investigation was mechanically identical to the NSTAR thruster. 4 A 
schematic of the ion thruster and internal discharge chamber wall probes is illustrated in Figure 1 . A screen 
grid and accelerator grid comprised the high-voltage ion extraction optics. The voltages at the screen and 
accelerator grid electrodes ranged from 650 V to 1500 V and -250 V to -180 V, respectively, as thruster 
input power varied from 0.5 kW to 4.6 kW. It should be noted that for this engine, the ratio of peak to 
average ion beam current density as measured in the plume varied from 0.45 to 0.55 over the range of 
conditions investigated. Additional details regarding the 30-cm engine, the power console, and the 
propellant feed system can be found elsewhere. 4,6,7 As can be seen in Figure 1, the aluminum discharge 
chamber was conical with a downstream cylindrical section. The magnetic circuit consisted of only three 
permanent magnet rings: 1) the cathode ring, located at the back-plate of the conical section; 2) the 
cylindrical section ring, located at the junction between the cylindrical section and the conical section; and 
3) the downstream ring, located at the most downstream flange of the cylindrical section. For these 
investigations, xenon propellant is used. 

The planar molybdenum wall probes used in this investigation each measured 6.4 mm in diameter. The 
wall probes were flush mounted with the surface of the anode at positions indicated in Figure 1 . The probes 
were designated according to the letters as indicated in Figure 1. Wall probe A was located near the 
midpoint between the cathode and cylinder magnetic cusps and wall probe B was located near the midpoint 
between the cylinder and pole piece magnetic cusps. The transverse magnetic field component across the 
surface of wall probes A and B was approximately 0.0021 T and 0.0023 T, respectively. The probes were 
electrically isolated from the anode wall and were held in place using modified compression fittings. The 
wall probes were biased relative to cathode potential using a variable DC voltage source to obtain 
Langmuir current-voltage characteristics from which plasma density and electron temperature were 
calculated. In order to obtain the average ion current density at the screen grid, the grid was negatively 
biased -25 V with respect to the discharge cathode. For this investigation, thruster discharge power ranged 
between 120 W and 465 W, which corresponded to a total thruster input power between 0.5 and 4.6 kW. 
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RESULTS AND DISCUSSION 


Background 

Between the cusps, the magnetic field lines run parallel to the surface of the anode. This field limits the 
cross-field diffusion of electrons to the anode surface. Figure 2a and 2b illustrate the variation in the 
normalized transverse magnetic field component at the anode as a function of position relative to the 
magnetic cusps. The normalized transverse field in the conical section of the discharge chamber is 
normalized to the transverse field intensity near the cathode cusp. For reference, the transverse magnetic 
field intensity at the anode one centimeter from the cathode ring magnetic cusp was 0.053 T. The 
normalized transverse field in the cylindrical section is normalized to the transverse field intensity near the 
cylindrical section magnetic cusp. The magnetic cusps are located at positions 0 and 18-cm of Figure 2a 
and positions 0 and 12-cm of Figure 2b. Along the anode, at distances greater than 2.5-cm away from the 
magnetic cusps, the transverse field does not vary appreciably. As can be seen in the figures, a significant 
portion of the anode surface lie in these regions of relatively constant transverse magnetic field. Indeed, 
over 70% of the anode's surface area lies in these regions of relatively constant transverse magnetic field. It 
therefore follows that the magnitude of the magnetic field at the anode between the cusps plays a critical 
role in determining discharge efficiency. The transverse magnetic field component in these regions is less 
than 20% of the transverse component measured near the cusps. Plasma collection at anode surfaces 
located close to the magnetic cusps is associated with magnetic cusp physics and is not the subject of this 
investigation. 

The focus of this work is to ascertain the effect that the inter-cusp magnetic field has on ion and 
electron transport to anode surfaces. Accordingly, this work is divided into two sections: 1) ion transport 
and 2) electron transport. 


Ion Transport 

Ion loss to the anode between the magnetic cusps and to the screen grid was determined by biasing 
wall probes A and B and the screen grid negative of the cathode potential to obtain the ion saturation 
current density. Figure 3 depicts a schematic of ion losses to bounding discharge chamber surfaces in the 
cylindrical section. The screen grid, which has a 67% open area, is the most downstream boundary of the 
discharge chamber’s cylindrical section. If preferential ion drift to the ion optics occurs in the NSTAR- 
derivative thruster, then it can be expected that the ion current density measured at the cylindrical section 
sidewall as measured by wall probe B should be less than the average ion current density measured at the 
biased screen grid. A plot of the ratio of ion current density measured at the screen grid to the ion current 
density measured at wall probe B is shown in Figure 4. For the discharge powers investigated, the ratio 
was nearly unity. This result is surprising in that it does not indicate a preferential flow of ions to the optics 
at any discharge power as has been suggested by an earlier study. 1 It should be pointed out that in the 
earlier study, the magnetic circuit as well as the anode geometry differed significantly from the NSTAR- 
derivative ion thruster studied in this present investigation. As determined from measurements presented 
here, the ion losses to the cylindrical section boundaries and the screen grid are approximately equal. It 
should be pointed out that such ion flow is consistent with the work of Masek. 8 

Past studies have suggested that the extracted ion beam fraction (defined as the ratio of the ion beam 
current to the sum of the ion current lost to discharge surfaces and the ion beam current) is considerably 
larger than that that would be expected from the ratio of screen grid physical open area to total discharge 
chamber surface area. 1,9,10 This result is attributable in part to an increase in screen grid transparency during 
high voltage extraction. 8,11 The remaining disparity can be addressed through an examination anode wall 
probe data. Figure 5 illustrates the behavior of the ion saturation current ratio of wall probe B to wall probe 
A as a function of discharge power. In all cases, the ion current measured at probe B was approximately a 
factor of 1.55-times larger than that measured at probe A. If the discharge ions escaping to the ion beam 
optics and to the cylindrical section walls obey the Bohm Criterion, the flow to the walls is at the ion 
acoustic velocity. Ion escape rates to anode surfaces should be proportional to the product of the square- 
root of the local electron temperature and the plasma density. Because the plasma density near the anode in 
both the conical and the cylindrical section of the discharge chamber was determined to be approximately 
equal, 10 the disparity between the ion wall flux measured at probe A and probe B must be due to the 
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electron temperature. In this regard, the ratio between the two probes' current densities should be 
approximately equal to the ratio of the square root of the electron temperature measured at the two probes. 
Plots of these ratios as a function of discharge power are shown in Figure 6. As can be seen from the plot, 
the average value of the ratio of the square root of the electron temperature is approximately 1.45, which 
within experimental error (10 to 15%) accounts for the difference in ion currents measured in the conical 
and cylindrical sections. This agreement supports the notion that a much higher percentage of ion 
production occurs just upstream of the ion optics and that the ion flow depends primarily on local plasma 
density and electron temperature. Because of the electron temperature disparity between the two discharge 
chamber sections, ion loss rates to the anode are largest in the cylindrical section, which in this case is just 
upstream of the optics. In this respect, ions are lost to the beam extraction optics at a rate larger than they 
are lost to anode surfaces in the conical section. Because ion losses to different sections of the discharge 
chamber are not uniform, the extracted ion fraction is not well represented by the ratio of screen grid open 
area to total discharge chamber surface area. 

The nature of ion flow suggests the importance of understanding plasma conditions, particularly the 
local electron temperature and density, within the discharge chamber. Care must be taken such that the 
discharge is not only localized near the optics, but also that the ratio of the open area of the screen grid to 
the anode surface area near the optics is sufficiently large. 


Electron Transport to Anode Surfaces between Magnetic Cusps 

The discharge chamber magnetic field controls the flow of electrons to the anode. Ideally, in a ring 
cusp discharge chamber, most of the electron current collection occurs at the magnetic cusps. Under these 
conditions, the magnetic field between the cusps minimizes the flow of energetic electrons to other anode 
surfaces. In practice, depending on the field strength between the cusps, a significant fraction of the 
discharge current can be collected at the anode between the cusps. It was found in an earlier NSTAR- 
derivative thruster study that up to 30% of the discharge current can be collected at anode surfaces between 
the cusps. 10 Because minimizing the flow of electrons to anode surfaces between the cusps improves 
thruster discharge efficiency, a understanding the nature of electron diffusion to anode surfaces between 
cusps is important. 

In general the radial flow of electrons to the anode in the presence of a magnetic field may be 
expressed as: 12 


J e =-e-fd-n e . E s -e-D ± 



( 1 ) 


Here mobility is defined as: 
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m • v e 


( 2 ) 


Of interest is the functional behavior of the cross-field diffusion coefficient. In general, plasma flow across 
a magnetic field is a complicated problem. Under some conditions, plasma flow obeys classical formulism 
in which the cross-field diffusion coefficient can be described by: 
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Here the electron diffusion coefficient under conditions where the magnetic field is zero is defined as: 


D 


e 


k-T. 

m ■ v e 


( 4 ) 


Under certain conditions, potential oscillations driven by plasma instabilities give rise to anomalous 
cross-field diffusion that is substantially larger than classical predictions. In a variety of cases, such 
anomalous diffusion is best described by the Bohm diffusion coefficient: 13-15 


D t = 


kT 


16-e-B 


(5) 


In this work, using Equation 1 and measured plasma properties at wall probes A and B, the electron 
current collected at and above plasma potential is calculated as a function of wall probe bias voltage. The 
total electron collision frequency, which is the sum of the electron-neutral collision frequency and the 
electron-ion collision frequency, was calculated based on ensemble averaged momentum exchange 
between the electrons and the background plasma using low energy electron-neutral and electron-ion cross 
section data. 16 Here the relaxation rate is determined by test particle (ion and neutrals) interactions with a 
Maxwellian electron distribution described by the measured electron temperature. 1 Neutral densities inside 
the discharge chamber were calculated from the propellant utilization efficiency. The calculation of the 
electron current is done assuming two separate cases: 1 ) Bohm diffusion and 2) Classical diffusion. The 
calculated electron current was then compared with the measured electron current to determine whether 
Bohm or classical best describes electron diffusion to anode surfaces between the magnetic cusps. 

In order to calculate the electron current at the wall probes as a function of probe potential, the electric 
field at the sheath edge as well as the electron density gradient must be determined. As illustrated in 
Figure 7, for all cases investigated in this work, the electron current at and above the apparent plasma 
potential was a linear function of the potential of the probe relative to the plasma potential. The constancy 
of the slope suggests the notion that the electric field in Equation 1 is proportional potential difference 
across the sheath. Equation 1 can then be recast: 


J e = -e • li • n e , ‘A V d eff - e-D 


dn e 

dr 


(D 


Here, the electric field at the sheath edge, may be expressed as the product of the sheath potential, AV , 
and the constant of proportionality, d cf j . The slope of the line described by equation 1\ 


— e-/U-?l e • d eff , can be obtained from a linear fit (see Fig. 7). From the slope, d e ^ can be calculated and 

thus be used to determine the effective electric field at the sheath edge. This effective electric field is a 
consequence of potentials that leak out of the sheath and give rise to electron drift toward the probe’s 
surface. Such fields cause the collected electron current to increase as a function of increasing sheath 
potential. 

The second term in equation (1') is associated with cross-field diffusion. Classically, in the presence of 
the transverse magnetic field, the characteristic diffusion step is a Larmor radius. 1 ' This scaling provides a 
simple means for estimating the electron density gradient. The cross-field spatial derivative of the electron 
density was taken to be the difference between plasma density in the bulk plasma (outside the sheath) and 


the density at the probe’s surface divided by the electron Larmor radius at the probe, 



\ 




Using equation (1'), the electron current was calculated for both the Bohm diffusion coefficient and the 
classical cross-field diffusion coefficient. Figures 8 and 9 illustrate the behavior of the ratio of the 
calculated electron current to the measured electron current as a function of sheath voltage for the two wall 
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probes at the 302 W discharge power condition (2.3 kW thruster input power) and at the 228 W discharge 
power condition (1.7 kW thruster input power). As can be seen from the figures, the electron collection 
between the cusps is best described by the classical cross-field diffusion coefficient. To within 
experimental error (approximately 35%), the classical relation predicts the magnitude of the electron 
current over a wide discharge power (200 to 465 W) range. Electron current calculated assuming a Bohm 
diffusion coefficient was over an order of magnitude larger than that which was measured. Indeed, for the 
magnetic field intensities present between cusps, if cross-field diffusion were Bohm-like, then the utility of 
magnetic containment would be defeated due to excessive electron flow to anode surfaces between the 
cusps. It should also be pointed out that the classical relation for the electron current also tracked the 
functional behavior of the electron current with increasing probe voltage. The rate of change in the electron 
current with increasing probe voltage as calculated using the Bohm diffusion coefficient case was less than 
that actually measured as evidenced by the monotonically decreasing ratios presented in Figures 8 and 9. 
Finally, it should be pointed out that the ratio of the electric field diffusion term to the magnetic field cross- 
diffusion term in equation 1 ranged from approximately 0. 1 to 1 for classical diffusion with increasing 
probe voltage. 

The fact that electron cross-field diffusion to inter-cusp anode surfaces is classical implies a number of 
interesting consequences: 1) Electron diffusion to the anode surface between cusps is severely reduced. 

2) Operation at higher propellant flow rates should increase electron losses to the anode because of the 
classical diffusion coefficient's dependence on the electron collision frequency. 3) Modest increases in 
magnetic field strength should lead to significant improvements in electron containment. 

Agreement between measured the electron current and the calculated classical current was best at the 
higher thruster power conditions (>1.4 kW). It should be pointed out that deviations from classical cross- 
field diffusion tended to increase at the lower power operating conditions. At these lower powers, the 
calculation for the electron current using the classical cross-field diffusion coefficient tended to under- 
estimate the measured value. The under-prediction trend at probe A ranged from 1% at the 2.3 kW input 
power condition to over 20% at the 1.4 kW input power. The discrepancy at wall probe B was somewhat 
larger than at probe A with the under-prediction trend ranging from 10% at the 2.3 kW input power 
condition to over 30% at the 1 .4 kW input power. 

Prior research suggests that deviations from the classical cross field diffusion may be due a plasma 
instability. 12 * 18 Additionally, ionization along the flux lines between the cusps may become more important 
with decreasing discharge power. Drainage of inter-cusp flux lines via cross-field diffusion or a plasma 
instability would also tend to increase the magnitude of cross field current collected at the probe. 10 These 
mechanisms may be necessary to sustain the magnitude of the discharge current at the lower thruster 
powers. The fact that measured electron diffusion to the anode between the magnetic cusps increases above 
the classical rate at the lower thruster power conditions suggests that electron containment between the 
cusps also degrades in concert. This degradation is consistent with the increased discharge losses at the 
lower power conditions. In this respect, characterization of electron diffusion phenomena to the anode is a 
key aspect for understanding discharge performance. 


CONCLUSIONS 

It was found that ion flow to anode surfaces just upstream of the ion optics is approximately equal to 
the average ion current density at the ion optics. Additionally, because ion losses to different sections of the 
discharge chamber are not uniform, the extracted ion fraction is not well represented by the ratio of screen 
grid open area to total discharge chamber surface area. Inter-cusp electron flow to the NSTAR-derivative 
ion thruster anode was found to follow classical diffusion predictions. This data suggests that modest 
increases in magnetic field intensity between the cusps can significantly improve electron containment. 
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Figure 1 . — Half-plane cross-section of NSTAR- 
derivative thruster revealing wall probes and 
magnetic circuit. 
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Normalized transverse magnetic field at anode Normalized transverse magnetic field at anode 



Location along conical section downstream of cathode cusp, cm 



Position along cylindrical section downstream of cylinder ring, cm 

Figure 2(a). — Normalized transverse magnetic field between the cathode ring magnetic 
cusp and the cylindrical section magnetic cusp. This profile illustrates the transverse 
field behavior in the conical section of the discharge chamber, (b) Normalized trans- 
verse magnetic field between the cylindrical section magnetic cusp and the down- 
stream magnetic cusp. This profile illustrates the transverse field behavior in the 
cylindrical section of the discharge chamber. 
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Figure 3. — Depiction of ion losses to the bounding discharge chamber anode 
surfaces and the ion optics. 



Figure 4. — Ratio of screen grid ion current density to ion current density measured at wall probe B 
as a thruster discharge power. 
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Figure 5. — Variations in the ratio ion saturation current measured at probe B to probe A as a function 
of discharge power. 
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Figure 6. — Variations in the square-root of the electron temperature ratio of wall probe B to wall 
probe A. 
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igure 7. — Electron current collected at wall probes A and B as a function of probe 
voltage above the plasma potential. 
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Figure 8. — Comparison between classical and Bohm electron flow to the anode 
between magnetic cusps at probe A (1 .7 and 2.3 kW input thruster power). 
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Figure 9. — Comparison between classical and Bohm electron flow to the anode 
between magnetic cusps at probe B (1 .7 and 2.3 kW input thruster power). 
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